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Abstract 

A non-stationary method for tunneling description of non-relativistic particles and 
photons through a barrier on the basis of consideration of the multiple internal reflections 
of vawe packets in relation of barrier boundaries is presented. The method is described in 
details and proved in the case of the one-dimentional tunneling of the particle through the 
rectangular barrier. For problems of the tunneling of the particle through the spherically 
symmetric barrier and of the photon through the one-dimensional barrier the amplitudes 
of transmitted and reflected wave packets in relation to the barrier, times of the tunneling 
and the reflection are found using of the method. Hartman's and Fletcher's effect is 
analysed. 
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1 Introduction 



The approach for description of a propagation of a nonrelativistic particle above a barrier using 
the account of multiple internal reflections of stationary plane waves between barrier boundaries 
which describe a motion of this particle in the region of the barrier, was considered in series of 
articles and is known for a long time PiJEJIHI- Thus, stationary solutions were studied only in 
the previous articles. 

To apply this approach for a solution of a problem, one can need in expressions for a wave 
function (w.f.) in the region of the barrier to separate components having fluxes, directed in 
opposite sides. For a problem of the particle propagating above the barrier it appears primely 
enough. So, considering an one- dimensional (ID) rectangular barrier, plane waves e ±thx can be 
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taken as such solutions, where k is a wave vector. To receive solutions using this approach for 
a problem of the tunneling of the particle under such barrier it appears more complicatedly, 
because in a consideration of the tunneling as a stationary process the decreasing and increasing 
components of the stationary w.f. in dependence on x (being analytic continuations of relevant 
expressions of the waves for the case of above-barrier energies) in a sub-barrier region correspond 
to zero fluxes separately and it is not correct to use them as the propagating waves from a 
physical point of view. If to define expressions for such waves by another way (for example, 
having required an existence of the nonzero fluxes in the barrier region at each step of this 
approach), then we receive a divergence of the expressions for the waves for the above-barrier 
and the sub-barrier cases. However, the flux calculated on a basis of a complete stationary 
w. f., is not equal to zero, and, therefore, the tunneling of the particle under the barrier exists. 

As a further development of a time analysis of tunneling processes submitted in articles 
|H El El here we represent the non- stationary solution method of a problem of tunneling of 
a nonrelativistic particle or a photon through a barrier using multiple internal reflections of 
fluxes in a sub-barrier region in relation to barrier boundaries (we name this approach as the 
method of multiple internal reflections). In given article we study an one- dimensional and a 
spherically symmetric problems. 

At analysing the tunneling (or the propagation) of the nonrelativistic particle, an important 
specific feature of this method is a description of a particle motion using non-stationary wave 
packets (w.p.). Due to this one can determine correctly the packets propagating in different 
directions in a barrier region, fulfil a time analysis of the particle tunneling (propagation) and 
study in details this process in an interesting time moment or in relation to a concrete point of 
space. For obtaining time parameters of the tunneling, this method has shown itself convenient 
and simply enough. 

For the problem with a spherically symmetric barrier the reflected and transmitted w.p. 
in relation to this barrier are propagate in one direction. Stationary solution methods do 
not allow to separate the w.p., transmitted through the barrier, from the w.p., reflected from 
the barrier. Using the method of multiple internal reflections, one can find amplitudes and 
expressions for these w.p.. In result, it appears possible to present S-matrix in a form of a sum 
of two components corresponding to stationary parts for the reflected and transmitted w.p. in 
relation to the barrier. The sum of the stationary parts for these w.p., obtained by this method, 
converges with the expression for a scattered wave obtained by an usual stationary method. 
The spherically symmetric problem with use of given approarch is considered for the first time. 

At first we consider the problem of the tunneling of the nonrelativistic particle through the 
one-dimensional rectangular barrier. This problem is a test one and allows to analyse specific 
features of this method. 

Further the problem of the tunneling of the particle through the spherically symmetric 
barrier, which radial part has a rectangular form, is solved. For it amplitudes of the transmitted 
and reflected w.p., total times of the tunneling and reflection in relation to the barrier are found. 
Hartman's and Fletcher's effect is analysed. An expression for S-matrix is presented in a form 
of a sum of two components corresponding to amplitudes of stationary parts for the transmitted 
and reflected w.p.. The time parameters using of the method of multiple internal reflections 
are found for the first time. 

One can apply the method to a problem solution of the tunneling of the particle through a 
spherically symmetric barrier, which radial part has an arbitrary form, if a general stationary 
solution for a w.f. is known for this potential. Some problems with various barrier forms are 
considered. The problems are selected so that to show better features of the method at their 
solution. 
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At a finishing of the article a possibility to use the method in a problem of a tunneling 
of photons through an one-dimensional rectangular barrier is studied. On a basis of a given 
analysis the method is proved for the problem with the photons. Using a found transforma- 
tion the results and solutions of the problem with the nonrelativistic particle transform into 
corresponding expressions for the problem with the photons. Hartman's and Fletcher's effect 
is analysed. 

2 Tunneling of a particle through an one-dimensional 
rectangular barrier 

Let's consider a problem of tunneling of a nonrelativistic particle in a positive x-dirrection 
through an one-dimensional rectangular potential barrier (see Fig. |IJ). Let's label a region I 
for x < 0, a region II for < x < a and a region III for x > a, accordingly. Let's study an 
evolution of its tunneling through the barrier. 

In a beginning we consider a standard approach to a solution of this problem pj Ej. Let's 
consider a case when levels of energy lay under than a height of the barrier: E < V%. 

The tunneling evolution of the particle can be described using a non-stationary consideration 
of a propagating w.p. 

+oo 

if>(x, t)= J g{E - E)^{k, x)e- iEt/h dE, (1) 
o 

where the stationary w.f. has a form: 

e %kx A R e~ tkx , for x < 0; 
ip(x) = { ae^ x + /3e~^ x , for < x < a; (2) 
A T e tkx , for x > a; 



and k = jrV2mE, £ = \^2m(Vi - E), E and m are the total energy and mass of the particle, 
accordingly. The weight amplitude g(E — E) can be written in a form of gaussian [3] and 
satisfies to a requirement of the normalization / \g(E — E)\ 2 dE = 1, value E is an average 
energy of the particle. One can calculate coefficients At, A^, a and (3 analytically, using a 
requirements of a continuity of w.f. tp(x) and its derivative on each boundary of the barrier. 

Substituting in Eq. ((TJ) instead of ip(k,x) the incident ip inc (k,x), transmitted ip tr (k,x) or 
reflected part of w.f. ip re f(k,x), defined by Eq. (J2j), we receive the incident, transmitted or 
reflected w.p., accordingly. 

We assume, that a time, for which the w.p. tunnels through the barrier, is enough small. 
So, the time necessary for a tunneling of an a-partile through a barrier of decay in a-decay of 
a nucleus, is about 10~ 21 sekonds 0. We consider, that one can neglect a spreading of the w.p. 
for this time. And a breadth of the w.p. appears essentially more narrow on a comparison with 
a barrier breadth jHEUHl- Considering only sub-barrier processes, we exclude a component of 
waves for above-barrier energies, having included the additional transformation 

g(E-E)^g(E-E)6(V 1 -E), (3) 

where ^-function satisfies to the requirement 



0{rj) 



for rj < 0; 
for rj > 0. 
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The method of multiple internal reflections considers the propagation process of the w.p. 
describing a motion of the particle, sequentially on steps of its penetration in relation to each 
boundary of the barrier Using this method, we find expressions for the transmitted 

and reflected w.p. in relation to the barrier. 

At the first step we consider the w.p. in the region I, which is incident upon the first (initial) 
boundary of the barrier. Let's assume, that this package transforms into the w.p., transmitted 
through this boundary and tunneling further in the region II, and into the w.p., reflected from 
the boundary and propagating back in the region I. Thus we consider, that the w.p., tunneling in 
the region II, is not reached the second (final) boundary of the barrier because of a terminating 
velocity of its propagation, and consequently at this step we consider only two regions I and 
II. Because of physical reasons to construct an expression for this packet, we consider, that its 
amplitude should decrease in a positive x-direction. We use only one item /3exp(— £x) in Eq. 
(J2J), throwing the second increasing item aexp(^x) (in an opposite case we break a requirement 
of a finiteness of the w.f. for an indefinitely wide barrier). In result, in the region II we obtain: 

+00 

tPI(x, t)= J g{E - E)6{V 1 - E)p°e^ x - lEt/n dE, for < x < a. (4) 


Thus the w.f. in the barrier region constructed by such way, is an analytic continuation of a 
relevant expression for the w.f., corresponding to a similar problem with above-barrier energies, 
where as a stationary expression we select the wave exp(ik 2 x), propagated to the right. 

Let's consider the first step further. One can write expressions for the incident and the 
reflected w.p. in relation to the first boundary as follows 

-00 



VwO, t) = J g(E - E)9(V 1 - E)e ikx - iEt / n dE, for x < 


+00 

J 




(5) 

ipl ef (x,t) = J g{E - E)9(V 1 - E)A° R e- ikx - iEt / h dE, for x < 0. 



A sum of these expressions represents the complete w.f. in the region I, which is dependent 
on a time. Let's require, that this w.f. and its derivative continuously transform into the w.f. 
(HJ) and its derivative at point x = (we assume, that the weight amplitude g{E — E) differs 
weakly at transmitting and reflecting of the w.p. in relation to the barrier boundaries). In 
result, we obtain two equations, in which one can pass from the time-dependent w.p. to the 
corresponding stationary w.f. and obtain the unknown coefficients (3° and A° R . 

At the second step we consider the w.p., tunneling in the region II and incident upon the 
second boundary of the barrier at point x = a. It transforms into the w.p., transmitted through 
this boundary and propagated in the region III, and into the w.p., reflected from the boundary 
and tunneled back in the region II. For a determination of these packets one can use Eq. Q 
with account (JHJ), where as the stationary w.f. we use: 

<fLc( k , x ) = <Plr(k,x) = pPe-* x , for < x < a, 

ip 2 tr {k,x) = A\e ikx , for x > a, (6) 

Lpl e j{k,x) = a°e^ x , for < x < a. 

Here, for forming an expression for the w.p. reflected from the boundary, we select an increasing 
part of the stationary solution a°exp(£x) only. Imposing a condition of continuity on the time- 
dependent w.f. and its derivative at point x = a, we obtain 2 new equations, from which we 
find the unknowns coefficients A% and a . 
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At the third step the w.p., tunneling in the region II, is incident upon the first boundary of 
the barrier. Then it transforms into the w.p., transmitted through this boundary and propa- 
gated further in the region I, and into the w.p., reflected from boundary and tunneled back in 
the region II. For a determination of these packets one can use Eq. (JTJ) with account Eq. (JHJ), 
where as the stationary w.f. we use: 

<Pinc( k > X ) = Vref{ k , X ), for < X < a, 

<pl(k,x) = A R e~ ikx , forx<0, (7) 
(Pr ef (k,x) = (3 l e~^ x , for < x < a. 



Using a conditions of continuity for the time-dependent w.f. and its derivative at point x = 0, 
we obtain the unknowns coefficients A l R and /3 1 . 

Analysing further possible processes of the transmission (and the reflection) of the w.p. 
through the boundaries of the barrier, we come to a deduction, that any of following steps can 
be reduced to one of 2 considered above. For the unknown coefficients a n , (3 n ,Aj, and A R , used 
in expressions for the w.p., forming in result of some internal reflections from the boundaries, 
one can obtain the recurrence relations: 

k + iC i£ + k i^ + k' ( a \ 

An on 2* , a . lka An+1 _ n 2ii W 



/|0 An on „—£a—ika An- 

A R ~ u i 77' T P 77 , 7 e ) A R 



a 



k + iC 1 r i£ + k ' ii + k 

Considering the propagation of the w.p. by such way, we obtain expressions for the w.f. on 
each region which can be written through series of multiple w.p.. Using Eq. (P) with account 
Eq. 0, we determine resultant expressions for the incident, transmitted and reflected w.p. in 
relation to the barrier, where one can need to use following expressions for the stationary w.f.: 

(p inc (k,x) = e ikx , forx<0, 

ip tr (k,x) = XM^e , for x > a, ^ 



n=0 
+oo 

(p ref (k,x) = £ A R e~ lkx , forx<0 

71=0 



Now we consider the w.p. formed in result of sequential n reflections from the boundaries of 
the barrier and incident upon one of these boundaries at point x = (i = 1) or at point x = a 
{i = 2). In result, this w.p. transforms into the w.p. ipl r {x,t), transmitted through boundary 
with number i, and into the w.p. il) l re j-{x,t), reflected from this boundary. For an independent 
on x parts of the stationary w.f. one can write: 

Vtr _m+ Vine Vref _ td+ Vine 



rp+ T me r rej _ j^-j- 



exp(— £x) exp{ikx) J exp{—ikx) exp{ikx)' 

2 2 2 2 

Vtr _ rp+ Vine Vref _ Vine {10) 

exp{ikx) 2 exp(— £x)' exp(£;c) 2 exp(— £r) ' 

V\r = T - Vine Vlef = R - V\ne 

exp{—ikx) 1 expire) ' exp(— £x) 1 exp{£x)' 

where the sign "+" (or "-") corresponds to the w.p., tunneling (or propagating) in a positive (or 
negative) x-direction and incident upon the boundary with number i. Using and Rf, one 
can precisely describe an arbitrary w.p. which has formed in result of n-multiple reflections, if 
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to know a "path" of its propagation along the barrier. Using the recurrence relations Eq. (jHJ. 
the coefficients and Rf can be obtained. 

An An+l 
T+ — RO T+ — T T~ — R 

J-\ -P, J-2 - T^T' J i - — ^T' 

P 11 



Using the recurrence relations, one can find series of coefficients a n , /3 n , A? and A R . How- 
ever, these series can be calculated easier, using coefficients Tf~ and Rf. Analysing all possible 
"paths" of the w.p. propagations along the barrier, we receive: 

i4k£e-t a - ika 



'-00 



E AVf = T+Tr 1+ E(RtRl 

n=0 V n=l 



+ OO / +OO \ knD 

ea r = Rt + T+R+Tr 1 + e (itf i*r) n = -|— 

n=0 V n=l / r s «6 



hex) 



2fc(i£ - k)e~^ a 



(12) 



where 



e « n = «° 1+ E(^r; 

n=0 V 71=1 

+ E> = /J »(i + + E( fl2+fll -)») = ?M|+* 

n=0 V i=l J r su b 



F sub = (k 2 -e)D„ + 2zk£D, 

D± = 1 ± e" 2 « a , 

, 9 , 9 „9 2m Vi 



(13) 



ft 2 ' 

All series J] a n , J] /5 n , J] and J] A R , obtained using the method of multiple internal re- 
flections, coincide with the corresponding coefficients a, (3, A T and A R of the Eq. (J2J), calculated 
by a stationary methods [UdE]]. Using the following substitution 

i£ - k 2 , (14) 



where k 2 = ^y2m(E — V\) is a wave number for a case of above-barrier energies, expression for 
the coefficients a n , (3 n , Aj, and A R for each step, expressions for the w.f. for each step, the total 
Eqs. (fT2*|) and (fH?|) transform into the corresponding expressions for a problem of the particle 
propagation above this barrier. At the transformation of the w.p. and the time-dependent w.f. 
one can need to change a sign of argument at ^-function. Besides the following property is 
fulfilled: 

+00 2 +00 2 

J2 a t + =i- (is) 



71=0 



R 

n=0 



3 Tunneling of the particle through a spherically sym- 
metric rectangular barrier 

3.1 Transmitted and reflected wave packets 

A problem of a motion of two interacting particles can be redused to a problem of one particle 
scattering in a spherically symmetric field. Let's assume, that the particle under an action of 
a central force 

— Vo, for r < Ri, (region I); 

V(r) = { V u for Ri < r < R 2} (region II); (16) 
0, for r > i? 2 , (region III). 
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is incident outside upon an external boundary of the barrier at point r = R2 (see Fig. 2). 

Let's study an evolution of tunneling of the particle through the barrier. We consider a case 
when the moment I = and levels of energy lay below than a barrier height. The tunneling 
evolution of the particle in time dependence can be described using a w.p. constructed on a 
basis of a stationary solution of the following form [Jj: 

^(r,M = — ( 17 ) 
r 

A(e~ iklT — e iklT ) , for r < R x , (region I), 

x(r) = { ae^ r + /?e~ 5r , for R x < r < R 2 , (region II), (18) 
e -ikr _|_ g e ikr^ f Qr r > (region III), 



where Yi m (9, <p) is a spherical function, k x = ^y2m(E + V ), £ = 2m(Vi — E), k = ^\/2mE. 
For the spherically symmetric problem in a case of sub-barrier energies we obtain: 

+00 

X (r, t)= J g(E - E)9(V l - E) X (k, r) e -' lEt / h dE, (19) 


V,(r) = V(r) + (20) 

where the second item in Eq. (f2T)|) is a centrifugal energy, which is equal to zero at I = 0, 
weight amplitude g(E — E) and average energy of the particle E are defined similarly to the 
one-dimensional problem (see Sec. |2J). 

At a stationary consideration of the solutions (fTHj) we describe the particle incident upon the 
external boundary of the barrier by a spherical wave exp(— ikr) convergent to the centre. And 
we describe the particle scattered on the barrier in the region III by a spherical wave S exp(ikr) 
divergent outside. The scattered wave takes into account both a possibility of a reflection of 
the particle from the barrier, which is written by the divergent wave, and a possibility of a 
penetration of the particle through the barrier, when in a beginning the particle tunnels from 
the region III to the region I, and then after some period of time it tunnels back from the 
region I to the region III and also is written by the divergent wave. Only one item S exp(ikr) 
contains the transmitted and reflected divergent waves, and it is impossible to separate them 
at the stationary consideration. 

As non-stationary, the method of multiple internal reflections allows to find a solution of this 
problem. Let's apply it to this problem. We study a propagation of a w.p. describing tunneling 
of the particle, sequentially on steps of its transmission in relation to each of boundaries of the 
barrier (similarly to the one-dimensional problem). In result of an analysis we come to a 
deduction, that any step in such viewing of the propagation of the w.p. along the barrier will 
be similar to one of 4 steps independent among themselves. Analysing these 4 steps further, 
one can obtain recurrence relations for finding coefficients A n , S n , a n and (3 n for an arbitrary 
step n. 

In result of multiple internal reflections (and transitions) in relation to the boundaries of 
the barrier a total time- dependent w.f. in each region can be written in a form of series, 
composed from convergent and divergent w.p.. Analysing possible "paths" of propagations of 
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these packets, one can calculate expressions for series of coefficients S n , A n , a n and (3 n : 
+00 1 

J2S n = T 2 -T 2 + (R^(l-RlR )+T l -R T l + ) = 

n=l ^sub 

u£ + h / 

+OQ T~T~ 4.,[h£ e -ik R 2+ikiRi-£(R2-Ri) 
j^n 1 2 S 

n=o F sub F sub (k + ig)(ki + i£) 

, 2k(l + b^ e 2ik 1 R l ) e -(Z+ik) R * ( 2i) 

, f - R{R V h + it, J y ' 



-00 







where 



(22) 



J2 a n = a 

n=o F suh F sub (k + i£) 

+00 

+~ _ So an ~ T2 n grq - RjRo) + t{RqT± 

n=0 -rt 2 P sub 

2k( — - ^ - e 2ifc i R i^ e 5(2Ri-«2)-jfeR2 
Vi^ + fci / 

/•;»/,(/•• • /X) 

= (f- J R+i? )(f-i? 2 + J Rr)-^r^i + = i+ , , 

_ (k - iQ(k 1 - iQ 2$(R 2 -Ri) _ k ~ ^ -2^(fi 2 -fli)+2ifcifli 

(fc + i^)(A;i + iO fc + i£ 

A; + if " ii + k 

T~ = — = C f£-Hfci)fli ij- = ^ = !i ~ fcl c 2 £^i 

a n + fci a n i£ + fci 

T - = i? = 1 (23) 

+ = /3 n+1 = 2fci (g+zfei)-Ri 73+ = = ^ ~ k l 2ik 1 R 1 

1 A" " i^ + ibi ' 1 A" i£ + £;i 

+ _ S^ 1 _ 2i£ {i+ik)R2 + _ a n+1 _ i£ - k 2 jR 2 

2 p" iZ + k ' 2 /3» + ^ 

where the coefficients and i?f are defined in relation to the boundary with number % {% — 
for r = 0, % = f for r = i?i and % = 2 for r = i^)- They can be calculated using of the 
recurrence relations between the coefficients S n , A n , a n and (3 n . 

Now we consider the incident, transmitted and reflected w.p. in relation to the barrier as a 
whole. Defining them for the region III, one can write: 

+00 _ 

Xinc(r, t) = J g(E - E)6(V 1 - E)e~ ikr - iEt l n dE, 


+00 _ 

Xtr (r, t) = J g(E - E)9(V 1 - E)S tr e^ Et ' h dE, (24) 


+00 _ 

X ref(r, t) = J g(E - E)9(Vi - E)S ref e^ Et ' h dE, 




where 



+00 

Str = J2 S n , 

Sref = S f (25) 
S = Sf r + S re f. 



The expression S represents a diagonal element of scattering matrix corresponding to the 
orbital moment I = 0. Thus, using the method of multiple internal reflections it appears 
possible to divide the S-matrix into two components corresponding to amplitudes of stationary 
parts of the transmitted and reflected w.p. in relation to the barrier as a whole. This property 
having physical sense, is obtained for the first time. 

The expressions for coefficients S n , A n , a n and (3 n for each step, the expression for the w.f. 
for each step, the coefficients and Rf, the series of the coefficients S n , A n , a n and (3 n under 
the substitution (JTljl (and also at replacement of a sign before argument for ^-function at a 
consideration of the non- stationary w.p.) transform into the corresponding expressions for a 
solution of a problem of a w.p. propagation above the barrier. The series (12 lj) of the coefficients 
S n , A n , a n and (3 n coincide with the corresponding coefficients S, A, a and (3 for Eq. (|TH|) . 
calculated by stationary methods. 



3.2 Tunneling and reflecting times in relation to the barrier 

One can determine an equation for a propagation of a maximum of the incident, transmitted 
and reflected w.p. in relation to the barrier for the spherically symmetric problem. For radial 
parts of non-stationary w.f. one can write: 

d d d 

— arg Xinc(r, t) = T^arg Xtr(r, t) = T^arg Xref(r, t) = const. (26) 

Let's consider the first step of the propagation of the w.p.. Let the w.p. is incident in 
the region III upon the external boundary of the barrier at point r = R 2 in a time moment 
t inc . Using Eq. we find the time moment t\ e j of leaving outside from this boundary the 
reflected w.p. in the region III: 

2mR 2 . <9arg S° . „ s 

^=«*.+-fir + *-5s- (27) 

Similarly, for a time moment t^ r of leaving outside from the external boundary of the barrier 
the n-multiple transmitted w.p. one can write: 

2mR 2 ^ <9arg S n , , 

c = w + _2 + a_s_. (28) 

Using Eq. (j2E|) at point r = R 2 , we find times necessary for the penetration of the total 
w.p. through the barrier (describing the tunneling of the particle through the barrier) and for 
the reflection of the w.p. from the barrier (describing the reflection of the particle from the 
barrier) : 



ph 2mR 2 darg S tr 

r tun itr tinc M +n QE , 

Ph + + 2mR i , ,dargS ref 



(29) 

—Ph / / " z _i_ t " 

ref "ref L inc j_ , t ii>~ 



hk dE 

For the problem of the w.p. tunneling under the barrier we receive 

Ph t d i^-h-iitl + hy 2 ^ 
t™ = h—aw 



^Ph 2m 

T ref 



dE & ^ + kf(ii + k 1 )F suh ' (30) 
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For the problem of the w.p. propagating above the barrier we write: 



Til = ti + h—ai. 



2m(R 2 -Ri) , d k 2 -h- (k 2 + h)e 2iklRl 



tun 



T ref — u ) 



hk 2 dE 6 (k + k 2 )(k 1 + k 2 )F above ' (31) 



where F above can be obtained from F sub using the substitution (JT3J). 

Let's consider a particle, which tunnels under a high enough and wide barrier. Then for the 
time of the tunneling we obtain the following expression (sequence of approaches: £,(R 2 — Ri) — > 
+oo, £ — > +oo, R 2 - i?i — ► +oo): 

Ph 2m 4mR 1 sin 2k 1 R i (l - 2 cos 2k 1 R l ) 

Ttun = m > + ^(1-^2^^) ' (32) 

The tunneling time does not depend on a width of the barrier (Hartman's and Fletcher's effect), 
but depends on k\ and R\. 

4 Tunneling of the particle through the spherically sym- 
metric barrier of a general view 

4.1 The particle propagates above the barrier 

In study of nuclear processes when a tunneling of partiles through a barrier is investigated, in the 
most cases the barriers of more complicated form than rectangular are used. So, the spherically 
symmetric two-humb potential of Strutinski has an enough important role in problems of fission 
and decai of nuclei. A degree of an exactitude of a description of the nuclear process depends 
on choice of a form of the potential. Therefore, we shall consider, as far as it is possible to 
use the method of multiple internal reflections for solving the spherically symmetrical problems 
with the barrier of a general view. 

Let's consider a particle propagating in a spherically symmetric potential field, which radial 
part has a barrier. Taking into account a behaviour of a radial part V(r) of the potential 
function in dependence on r, we divide the area of its definition r 6 [0; +oo[ on n regions. In 
each region let's replace the potential function V(r) by a function most close describing V(r) 
and for which an exact solution of the stationary Schrodinger equation exists (see Fig. ??). 
Passing to the problem of the particle propagation in the field of these approximated potential 
functions, we write the general solution for stationary w.f. in the form ()17)) . where its radial 
part can be written as 

{Aiai(k,r) + Bibi(k,r), for < r < R\, (region I), 
Ai(ii(k, r) + Bibi(k, r), for R^ < r < R u (region i), (33) 
A n a n {k } r) + B n b n (k, r), for r > i? n -i, (region N), 

where k = y i \ / 2mE, <Xi(k,r) and bi(k,r) are the partial solutions of the radial part of w.f. in 
region i, Ai and Bi are the normalization constants. 

Let's find the transmission and reflection coefficients of particle in relation to the barrier, 
and also the times necessary for transmission and for reflection of the particle in relation to the 
barrier, using the method of multiple internal reflections. To apply the method to this problem, 
one can need to present the general stationary solution of w.f. in each region in the sum of 
divergent and convergent vawes. 
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Using the Fourier transformation, one can write: 

Oi(k,r) = a~(k,r) + af(k,r), 
bi(k,r) = bt(k,r) + bf(k,r), 



(34) 



where 



o R l 

a~{k,r) = -= f dq f a^r^e^-^dr' , 



2,7V 



hoo Ri 



af(k,r) = ^= f dq f a(k,r')e tq{r - r,) dr' , 

\/27r J J 

lU-i 

Ri 

b~(k,r) = -= j dq ( biikyy^-^dr', 



2tT 



(35) 



-OO 
hoo Ri 



bf{k,r) = f dq f bifar'y^-^dr'. 

\/27r J J 


Taking into account the transformation 

Ci(k,r) = a~(k,r) + -^bt(k,r), 

A . Al 
4( k , r ) = -£at(k,r) + bf(k,r) 

one can write the general solusion (|3~3|) as 



(36) 



AxCi(k, r) + BiCi(k, r), for < r < R±, (region I), 
X(r) = { AiC~{k } r) + Bicf(k,r), for < r < R h (region i), (37) 
A n c~(k,r) + B n c+(k,r), for r > R n -\, (region N), 

In result, the general solution in every region i is represented as the sum of convergent 
vawes c~(r) and divergent vawes cf(r) (so, in case of a rectangular barrier in the region i such 
expressions equal to e~ lkiT and e iklT , accordingly). On the basis of these expressions using Eq. 
(fE?j) one can construct the non-stationary convergent and divergent w.p.. Writing the genegal 
solution for w.f. in every region in the form of linear combination of convergent and divergent 
w.p., one can apply the method of multiple internal reflections for solving the problem. 

At first we study the case, when the general stationary solution for w.f. in every region can 
be written uniquely as sum of convergent and divergent vawes cf (r). Then using of the method 
of multiple internal reflections for solving the problem, one can find the incident, transmitted 
and reflected w.p. in relation to the barrier at whole, and total w.p. in every region. It is 
enough convenient to use the coefficients Tf 1 and Rf (as in the spherically symmetric problem 
with rectangular barrier). We define these coefficients in relation to the boundary with the 
number i by such way (for the step j): 

B j i+l = T?B{; At 1 = RfBl ^ 
One can calculate these coefficients at consideration of first 2n + 1 steps: 

To - = 0; * = -gg = fi, (39) 
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T~ 

1 i 



TP 



dc i+1 (r] 
dr 



dcj +1 (r) 
dr 



dr 
dc i+1 (r] 

dr 



c A r 



dr 
(r) 



--it, 



dr 



dd 



dr 
dr 



-c- +1 {r)- Ci 



-c, r - c r 



dc? +1 (r) 

dr 
d c i+1 (r) 

dr 



dcf +1 {r) 

dr 
dcf(r) 



dc; (r) 



C; (r 



--Ri 



(40) 



dr 



-cj+i\ r ) 



4(i 



dr 
dc i+1 (r) 

dr 



dcf +1 (r) 



4+i( r 



d^ (r) 



-Ri 



dr ' ' dr 

Further using the method of multiple internal reflections, one can calculate the incident, 
transmitted and reflected w.p. in relation to the barrier. On the basis of these w.p. one can 
find the transmission and reflection coefficients and also the transmission and reflection times 
in relation to the barrier. Thus the transmitted and reflected w.p. can be written through St r 
and S re f, which sum is the diagonal element of the scattering matrix S at the orbital moment 
I. 



The value A n can be obtained from the normalization condition: 



-1/2 



\c n (k,r)\ 2 dr 



(41) 



Now we study the case, when the partial solutions of w.f. in some regions are not the 
convergent and divergent vawes. In representations and ()37j) one can need to know the 
values Bi/Ai. In this case at solving the problem we consider first 2n + 1 steps. Let the general 
solution for w.f. in first region be expressed through ai and b\. Analyzing the reflection of w.p. 
from point r = 0, one can obtain: 



0; i?o 



ai(0) 

Mo) 



Bi 



(42) 



If the w.f. in the first region is determined through cf uniquely, then it is need to use Eq. 
instead of Eq. (j4*2^) . Calculating the value Bi/Ai, then one can find the functions cf . 
Using the continuity condition for w.f. and its derivative in all boundaries between regions, one 
can find the recurrent relation for values Bi/Af. 



B 



Mr) 



daj 



dr 



dfjjr ) 
dr 



-Ci; 



Ai 



i+l 



i+l[ r ) 



Mr) 



" dMr) 
dr 

ai{r) + bi(r)^ 



Mr) 



db i+1 (r) 
dr 



-Ri 



(43) 



-Ri 



Having the values Bi/Ai, we obtain the convergent and divergent vawes cf (r) in every 
region. Then the solution of problem is fulfilled as in the previous case. 
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Applying the approach considered above for the solution of problem of particle propagation, 
when the potential V{r) is defined only in two regions (n = 2), one can find the incident, 
transmitted and reflected w.p. in relation to the barrier using Eq. (JT§j) for above-barrier region, 
where the radial parts from corresponding stationary w.f. have the form (at r > R±) 

X mc (r) = A 2 c 2 (k,r), 

Xtr(r) = S tr ct(k,r) = A j lT ^ + R " ct(k,r), (44) 

1 — ill 

Xref{r) = S ref ct(k,r) = A 2 R±c$(k,r). 

We find the coefficients Tf 1 and Rf from Eq. (|4U|). Using Eq. ()26|) for external boundary, 
one can obtain the times necessary for transmission and for reflection of particle in relation to 
the barrier. In result, we receive: 



Ttun = (arg B 2 c t( k , R l) ) ~ arg A 2C 2 (k, Ri) 



T, 



ref 



d T x T x R 

h— — arg x h At: 

dE & 1 - RtR ' 

^^( arg (^24(^,^1)) - ar S A 2 c 2 (k,R 1 



(45) 



d 

= h -Q^ sx ^> R i + Ar > 

At = h dE^^iR7y 

For the problem solution when potential V(r) is defined on three regions (n = 3), the 
expressions for radial parts of stationary w.f., describing the incident, transmitted and reflected 
w.p. in relation to the barrier, look like (at r > R 2 ) 

Xincir) = A^(k,r), 

(\ q \ a r 2 -r+(i?r(i - RjR ) + t 1 -r t 1 + ) 

Yf r (r) = b tr ci(k,r) = A?- t w t ; ttz T-rcT(k,r) (4o) 

Xref(r) = S ref ct(k,r) = A 3 R 2 c£(k,r). 

The transmission and reflection times of particle in relation to the barrier has the form 
(they are calculated at r = R 2 ) 

_ d T 2 -T 2 + (i?r(l - RfRo) + TT^Tf) 

dE g (1 - RfR )(l - RtRi) - R 2 T{ RqTi + ' 
d 

Tref = ^^g ar S R 2 + Ar ! ( 47 ) 

+ d ct(R 2 ) 



As an examples of method application we consider two problems. 
The particle propagates above the barrier of form (see Fig. ??) 

!— Vo, for < r < Ri, (region I); 
7 r D , • TT\ ( 48 ) 

— , tor r > Hi, (region 11). v ' 

r 
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We consider the case / 7^ 0. One can obtain the incident, transmitted and reflected w.p. in 
relation to the barrier from Eqs. (jUJ) and (|19|h taking into account the sign before argument of 
^-function for above-barrier energies, and transmission and reflection times from Eq. At 
consideration the first three steps of w.p. propagation along the barrier we find the coefficients 
T? and Rf using Eqs. (|4f)|) for n = 2. In the solutions one can need to fulfil the substitution 

cl(k,r) = Xkxi(r), &[(k,r) = Gfa p) - iFfa, p), ^ 

where 



k x = y2m(E + V ), 
pvk 

p = k(r), 

Hi ' (r) is the function of Hankel of 1-st and 2 sort, 

Gi(r],p) and Fi(r),p) are the irregular and regular Coulomb functions jTTj. The normalization 
constant A 2 can be obtained from Eq. (|4ip. 

Now we consider another problem when the particle propagates above the barrier of following 
form: 

{ar 2 — Vo, for < r < R±, (region I); 
7 t D / • T n (51) 

— , tor r > iii, (region 11). v y 

r 

Let's study the case 1 = 0. In the beginning we consider region I. The partial solutions 
for the radial part of stationary w.f. are the parabolic cylinder functions [llj: D u (±gr) and 
D u (±igr), where g = (8ap/h ) 1//4 . For the description of above-barrier motion of particle we 
choose the first two solutions D u (±gr), which are independent if v is non-integer. (Note that 
one can use the Whitteker's functions as such two independent solutions JT]. But these two 
functions can be presented in the form of linear combination of the parabolic cylinder functions 
D u (±gr).) Each of partial solutions can be presented in the form of sum of convergent and 
divergent waves: 

D v {±gr) = D;(±gr) + D+(±gr), 



u -n.1 

DJ±gr)~ = -= [ dq [ D v (±gr')e tq{r - r ^dr\ 

x/27T J J (52) 



-00 

-00 i?i 



D+(±gr) = —= f dq f D v (±gr')e iq{r - r,) dr' 
\/2n J J 



2tt 



Using such w.f., one can apply the method of multiple internal reflections to the solution 
of problem. In result, we find the incident, transmitted and reflected w.p. in relation to the 
barrier from Eqs. (}44"j) and (|19j) . taking into account the sign before argument of function for 
above-barrier energies, and transmission and reflection times from Eq. (|45jl . The coefficients 
T/* 1 and Rf can be obtained from Eqs. (|4()jl and (|42j) for n = 2 at substitution 

ai(fc,r) = D u (k,r), c£(k,r) = G (r], p) - iF (r], p), , . 

bxi^r) = D v (-k,r), c+(k,r) = G (v, p) + iF ( V , p), K } 

where r\ and p are defined in Eq. (jHUj) . Go(r],p) and F Q (r),p) are the irregular and regular 
Coulomb functions at I = 0. 
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4.2 The particle tunnels under the barrier 



Now we consider the problem of tunneling of particle under the barrier of spherically symmetric 
potential field. And the radial part of this barrier has a general view (see Fig. ??). 

Dividing the range of definition r e [0; +00 [ for potential V(r) on n regions, on each of them 
we approximate V(r) by function most close to it, for which there are the general solutions of 
w.f. for stationary Schrodinger equation. We divide the whole range so that the processes of 
sub-barrier tunneling and above-barrier propagation laid in the different regions. 

For regions, in which the energy levels considered by us, lay above the potential function 
V(r) (the particle propagates above the potential), the stationary solution for w.f. is represented 
as Eq. (J37)) (if necessary using the transformations (JMJ), (jBSj) and (JHSJ))- 

For regions, in which the viewed energy levels lay under the potential function V(r) (the 
particle tunnels under the potential), in the beginning we find the general solution for stationary 
w.f., assuming, that the energy levels lay above the potential function. One can need to present 
the general solution for w.f. as Eq. (pT7jl . separate the components corresponding to fluxes, 
directed to the opposite sides. Everywhere in expressions for w.f., where the property 



is used, one can need to redefine this expression for E < Vi(r), having changed the sign. So, 
in case of constant potential in dependence on r we obtain the Eq. (JHJ). Such substitution 
gives the following property: the resultant expressions for w.p. and also for stationary and 
non-stationary w.f. for the problem of tunneling of a particle under the barrier are the analytic 
continuation of the relevant expressions for a similar problem, when the particle propagates 
above the barrier. 

Having defined the expressions for stationary w.f. by such way, one can construct the 
relevant them w.p. on each region and apply the method of multiple internal reflections to 
solution of the problem. The further approach for obtaining the resultant expressions for 
incident, transmitted and reflected w.p. in relation to the barrier and also the times of tunneling 
and reflection differs by nothing from the approach for the problem solution in the above-barrier 
case. 

As an example, we consider the problem of tunneling of particle under the barrier (|48|) (see 
Fig. ??). We consider the case I 7^ 0. We divide the region II on two at point r = R 2 , which 
defines by requirement E = Vi{Rz). One can find the incident, transmitted and reflected w.p. 
in relation to the barrier from Eqs. (fTTH) and (|4fij) . and the times of tunneling and reflection 
from Eq. (|47jl . Analysing the first 5 steps of w.p. propagation along the barrier, we find the 
coefficients T i and Rf using the Eqs. (|4*0|) for n = 3. In these expressions on can need to fulfil 
the substitution 




(54) 




(55) 



where k, ki, r), p, and also Xkii( r )i Gi(ViP) an d Fi(v->P) are defined earlier. 
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5 Evolution of photon tunneling through one- 
dimensional undersized rectangular waveguide 

We use the analogy between photon and particle ID propagation and tunneling which consists 
not only in the formal mathematical analogy between the solutions of the time-dependent 
Schrodinger equation for nonrelativistic particles and of the time-dependent Helmholtz equation 
for electromagnetic waves but also in the similarity of the probabilistic interpretation of the 
wave function for a particle and of a an electromagnetic wave packet being the wave function 
for a single photon [Sj. For a hollow rectangular waveguide with variable section (like that used 
in the Cologne experiment ^2], see Fig. ??). The time-dependent wave equation for A, E, H 
(A is the vector potential with the subsidiary gauge condition div A = 0, E = —(l/c)dA/dt 
is the electric field strength, H = rot A is the magnetic field strength) is 

1 <9 2 A 

For boundary conditions (see, for instance, [3]) 

E y = for z = and z — a, , _^ 

E z = for y = and y = b 

the solution of the Eq. (|57|) can be represented as a superposition of the following monochro- 
matic waves: 

E x = 0, 

Ey = Eq sin (k z z) cos (k y y) exp [i(wt ± 7a;)] , (58) 
Ef = — E Q (k y /k z ) cos ik z z) sin (k y y) exp [i(wt ± jx)], 

where k 2 + k 2 +'-f 2 = w 2 /c 2 = (2ir/\) 2 , k z = ran /a, k y = nir/b, m and n are the integer numbers 
(for definiteness we have chosen the TE- waves). Thus, 

7 = 2tt[(1/A) 2 - (1/A C ) 2 ] 1/2 , (1/Ac) 2 = (m/2a) 2 + (n/2b) 2 , (59) 

where 7 is real (7 = Rej) if A < A c and 7 is imaginary (7 = i£ em ) if A > A c . Similar expressions 
for 7 were obtained for TH- waves [S]. 

Generally the non-stationary solution of Eq. (J56j) can be written as a wave packet constructed 
on the basis of monochromatic solutions (p)%j) . similarly to the solution of the time-dependent 
Schrodinger equation for nonrelativistic particles in the form of a wave packet constructed 
from monochromatic terms (for the problem of particle propagating above the ID rectangular 
barrier). Moreover, in the representation of primary quantization the probabilistic single-photon 
wave function is usually described by a wave packet (for instance, see E] and the relevant 
references therein) like 

A(r,t)= / ^K(k)A(k,r)e- iM , (60) 



j k 

k >0 

where A(\t, r) = exp (zkr) for propagation in vacuum and A(k, r) = (p(x) exp(ik y y +ik z z) with 

e ik x x _|_ aRe -ikxx^ region I, 
(p(x) = { ae~^ emX + (3e^ emX , region II, (61) 
a T e tkxX , region III 

2 

for propagation in the waveguide (Fig. ??). Here, r = (x,y,z), K(k) = J2 -^i(k)ej(k), e^e^ = 

i=l 

6ij, ej(k)k = 0, i,j = 1,2 (or y, z if kr = k x x), k = w/c = e/hc, |k| = k , Ki(k) is the 
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amplitude for the photon with impulse k and polarization i, and |f^(k)| 2 dk is then proportional 
to the probability that the photon has the impulse between k and k + dk in the polarization 
state ej. 

Though it is not possible to localize photon in the direction of its polarization, nevertheless, 
in a certain sense, for the one- dimensional propagation it is possible to use the space-time 
probabilistic interpretation of Eq. (ffiUJ) along axis x (the propagation direction) . It can be 
realized from the following. Usually one uses not the probability density and probability flux 
density with the corresponding continuity equation directly but the energy density s and the 
energy flux density s x (although in general they represent components of not a 4-dimensional 
vector but the energy-momentum tensor) with the corresponding continuity equation jH] which 
we write in the two-dimensional (spatially one-dimensional) form: 

< 62 > 

where 

EE + HH oRefEHL ,„* 
so = 5 , s x = (pi) 

87T Z7T 

and axis x is directed along the motion direction (the mean impulse) of the wave packet ([60)1 . 
Note, that for the spatially one-dimensional propagation the energy-momentum tensor of the 
electromagnetic field reduces to the two-component quantity — to the scalar term sq and 1- 
dimensional vector term s x for which continuity equation (}62|) is Lorentz-invariant. Then, as 
a normalization condition one chooses the equality of the spatial integrals of so and s x to the 
mean photon energy and the mean photon impulse respectively or simply the unit energy flux 
density s x . With this, we can define conventionally the probability density 

Pemdx = , 5*0 = / s dydz (64) 



J S dx' 

for the photon to be found (localized) in the spatial interval (x, x + dx) along axis x at the 
moment t, and the flux probability 

J e m,xdt = , S x = J s x dydz, (65) 

for the photon to propagate through point (plane) x in the time interval (t, t+dt), quite similarly 
to the probabilistic quantities for particles. Hence, in a certain sense, for time analysis along 
the motion direction, the wave packet ()60|) is quite similar to a wave packet for nonrelativistic 
particles and similarly to the conventional nonrelativistic quantum mechanics, one can define 
the same form of time operator as for particles in nonrelativistic quantum mechanics and 
hence the mean time and the distribution variance of times of photon (electromagnetic wave 
packet) passing through point x in both time and energy representations) jHj- Then, the 
same interpretation one can use for the propagation of electromagnetic wave packets (photons) 
in media and waveguides when reflections and tunneling can take place — in particular, for 
waveguides like depicted in Fig. ?? with spatially decreasing and increasing waves in Eq. (JBTJ. 
The only difference is in the momentum-energy relation (quadratic for particles and linear for 
photons). 

So, from rather simple calculations of J em ,x using Eqs. (jHOJ) - (JHSJ) , and using the given above 
definitions of E and H (see also [6|), one can obtain the following relation: 

J em , x = ReF(y,z)^p(x)^^^j, (66) 
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where the function F(y, z) depends on the boundary conditions of the waveguide (see Fig. ??) 
and calculated in jB] . Therefore under boundary conditions the flux density for photons can be 
obtained from the flux density for particles by simple replacing (—ih/2m) by F(y, z)[— iw/Ait]. 
At this substitution all results and relevant expressions (approach to the solution of a problem 
on the basis of consideration of multiple internal reflections of fluxes in the region of the barrier, 
phase tunneling and reflection times and other results), obtained above for the description of 
tunneling evolution of the particle through the barrier, also take place at the description of 
photon propagation. 

In the particular case of quasimonochromatic wave packets, under the same boundary con- 
ditions as considered for the problem of tunneling of a particle through ID rectangular barrier, 
we obtain the identical expression for the phase tunneling time: 

<n,em = j- for £ em L » 1. (67) 

From Eq. (jfi7j) one can see that when £ em L > 2 the effective tunneling velocity 

v tll = —ph (68) 

' tun,em 

is more than c, i.e. superluminal. This result agrees with the results of the microwave-tunneling 
measurements presented in |12j . 

Note, that for sub-barrier energies the nonlocality of a barrier as a whole takes place not 
only for nonrelativistic particles but also for photons. This property is the physical cause of 
the superluminality during the tunneling. 



6 Conclusions 

In this work the method of multiple internal reflections describing the process of tunneling of 
a nonrelativistic particles and photons through barriers of the various forms is presented. This 
method is the further development of a series of articles jHEJE], devoted to the time description 
of tunneling through a barrier. It uses the essentially non-stationary approach constructed on 
the basis of multiple reflections (and transmissions) of w.p. in relation to the boundaries of 
barrier. In result one can describe in dependence on time the process of tunneling of total 
w.p. describing the considered nonrelativistic particle or photon, through barrier and to study 
specific features of process in any interesting moment of time or in any point of space in details. 

The possibility of time description of tunneling through a barrier is one of the principle 
perspectives of this method in comparison with stationary approaches. 

The stationary one-dimensional problem of tunneling (and propagation) of a nonrelativistic 
particle through a rectangular barrier with accounting of the multiple internal reflections was 
earlier solved P3 121 El • For sub-barrier energies the plane waves in the barrier region (on the 
basis of which the complete expressions for w.f. were found) had zero fluxes. According to 
the physical understanding there is a problem of applicability of such approach to the problem 
solution. In the given article the substantiation of this approach is given on the basis of using 
the non-stationary w.p.. For this problem (being the test one) the phase time of tunneling 
and reflection in relation to the barrier at whole under solving the problem on the basis of the 
method of multiple internal reflections are introduced. 

Using the method of multiple internal reflections the problem of tunneling of a nonrelativistic 
particle through a spherically symmetric barrier is solved for the first time. Here, using this 
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method it is possible (as against the known stationary approaches) to separate the wave packet, 
transmitted through the barrier and describing a particle after its leaving outside after double 
tunneling through barrier, from the wave packet, reflected from the barrier describing a reflected 
particle (both packets are spherically divergent). For the diagonal element of scattering matrix 
with orbital moment / the following property 

qI ol I cl 

is fulfilled, i.e. the S'-matrix consists of two components corresponding to the transmitted and 
reflected wave packets in relation to barrier. This property has physical sense and is proved 
mathematically. 

We suppose, that the method will allow to describe such properties of nuclear processes, 
which are not explained by stationary methods. So, some experiments performed recently, 
have caused an increased interest to a bremsstrahlung in an a-decay of heavy nuclei . This 
phenomenon is interesting in that includes both a radiation of photons in a propagation of an 
a-partile in an electromagnetic field of an daughter nucleus, and a tunneling of the a-partile 
through the decay barrie. Now the effect of the photon radiation in the tunneling of the a- 
partile under the barrier is investigated unsatisfactoryli. For a description of this process some 
stationary methods allowing to calculate a spectrum of the bremsstrahlung are created. But 
in comparison with the experimental data one can see that each stationary approach describes 
the phenomenon with a small degree of an exactitude. Besides the minima and maximas are 
registered in the spectrum for some nuclei, while the stationary methods give a monotonically 
decreasing curve for the spectrum. We assume, that based on a space-time approach the method 
of multiple internal reflections will allow to explain the peaks in this spectrum. A preliminary 
analysis shows, that these peaks correspond to resonance levels of the a-decay of the researched 
nucleus and they can be evaluated using of the method. 

In article the possibility to apply the method for ID problem of photon tunneling through a 
rectangular barrier is explored. On the basis of the given analysis the analogy (having a math- 
ematical substantiation and physical sense) between wave packets (and also between problem 
setting, boundary conditions) describing both propagation and tunneling of a nonrelativistic 
particle and photon, is shown. In result, it is possible to apply the method of multiple internal 
reflections for the problem with photons for the first time. At the found transformation the ob- 
tained results for the problem of particle tunneling through a barrier transform into the relevant 
expressions for the problem of tunneling of photons. The tunneling durations are found. For 
enough wide (and high) barrier there is an effect of propagation of wave packet with velocity 
more, than velocity of light (Hartman's and Fletcher's effect). 

The superluminal phenomena, observed in the experiments presented in [12] and later in 
other papers (for example, see the relevant references in El EB] ) , generated a lot of discussions 
on relativistic causality. And in connection with this, also an interest for similar phenomena, 
observed for the electromagnetic pulse propagation in a dispersive medium |14j . was revived. 
The known way of usual understanding consists in explaining the superluminal phenomena 
during tunneling on the base of a pulse attenuated reshaping (or reconstructing) discussed at 
the classical limit earlier by GSj: the later parts of an input pulse are preferentially 

attenuated in such a way that the output peak appears shifted toward earlier times, arising 
from the forward tail of the incident pulse in a strictly causal manner [T7] . 

Already for long there was ascertained that the wavefront velocity of the electromagnetic 
pulse propagation, when pulses have a step-function envelope, cannot exceed the velocity of 
light c in vacuum [T5| ITH] . Namely in this the principal demand of the relativistic (Einstein) 
causality consists . This conclusion was confirmed by various methods and in various processes, 
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including tunneling ^Hl E3 EDI EH 1221 123 121]- Note, that it is known from the momentum- 
energy Fourier- analysis of an electromagnetic wave packet with the step-function form of the 
forward edge, that such a wave packet contains components with large (up to the infinite) 
energies, i. e. above-barrier energies, for which the superluminality is absent. 

One of the argued now problems consists in the absence of a step-function form of forward 
edges for realistic wave packets [TBI EH1 121] • In such cases the conclusions of [JJ3 HH] can seem 
to be inapplicable. An interesting approach to analyse the form of causality namely in such 
cases was proposed in f2~o) . 

Finally, from the analysis of first step in solving the problem by method of multiple internal 
reflections one can see that the tunneling process at sub-barrier energies is a non-local phe- 
nomenon becauce during tunneling the entering w.p. fills up the whole barrier atonce and w.p. 
feels immediately both barrier walls (boundaries). 
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Figure 1: One-dimensional rectangular barrier 
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Figure 2: Spherically symmetric rectangular barrier 
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